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Abstract Following an extremely interesting idea [I], published long ago, the work function associated 
with the emission of ultra-relativistic electrons from magnetically deformed metallic crystal of astro- 
physical relevance is obtained using relativistic version of Thomas-Fermi type model. In the present 
scenario, surprisingly, the work function becomes anisotropic; the longitudinal part is an increasing 
function of magnetic field strength, whereas the transverse part diverges. 
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1 Introduction 

In condensed matter physics the work function is the minimum energy needed to remove an electron 
from a solid to a point immediately outside the solid surface or equivalently, the energy needed to move 
an electron from the Fermi level into vacuum. Here the word immediately means that the final electron 
position is far from the surface on the atomic scale but still close to the solid on the macroscopic scale. 
It is also known that the work function is a characteristic property for any solid face of a substance 
with a conduction band, which may be empty or partly filled. For a metal, the Fermi level is inside the 
conduction band, indicating that the band is partly filled. For an insulator, however, the Fermi level 
lies within the band gap, indicating an empty conduction band; in this case, the minimum energy to 
remove an electron is about the sum of half the band gap, and the work function. 

In the free electron model the valence electrons move freely inside the metal but find a confining 
potential step, say C at the boundary of the metal. In the system's ground state, energy levels below 
the Fermi energy are occupied, and those above the Fermi Level are empty. The energy required to 
liberate an electron in the Fermi Level is the work function and is given by Wf — C — fx e , where fi e 
is the Fermi energy. Therefore, the work function of a metal is usually defined as the smallest energy 
needed to extract an electron at zero temperature. Formally, this definition is made for an infinitely 
large crystal plane in which one takes an electron from infinitely deep inside the crystal and brings it 
through the surface, infinitely far away into the vacuum. The exact definition of work function, which is 
valid in the atomic scale is the work done in bringing an electron from far below the surface, compared 
to atomic dimensions but not far compared to crystal dimensions. 

The work function of metals varies from one crystal plane to another and also varies slightly with 
temperature. For a metal, the work function has a simple interpretation. At absolute zero, the energy 
of the most energetic electrons in a metal is referred to as the Fermi energy; the work function of a 
metal is then equal to the energy required to raise an electron with the Fermi energy to the energy 
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level corresponding to an electron at rest in vacuum. The work function of a semiconductor or an 
insulator, however, has the same interpretation, but in these materials the Fermi level is in general 
not occupied by electrons and thus has a more abstract meaning (see [2J|3J|IJ[5J|1>1[7] as some of the 
fundamental papers on work function). 

The work function is associated with three types of electron emission processes: photo-emission, 
thermionic emission and field emission. All of them have a large number of applications in various 
branches of science and technology; starting from electronic valves, a very old type electronic device, 
to modern opto-electronic devices used to convert optical signals to electrical signals etc. The other 
important applications are in photo- multiplier tube (PMT), CCD, etc. In the recent years, it is found 
that the field emission process has a lot of important applications in modern nano technology. 

Similar to the down to earth application of work function in condensed matter physics and its various 
scientific and technological applications, the work function also plays vital role in many astrophysical 
processes, e.g., in the formation of magneto-sphere by the emission of electrons from the polar region 
induced by strong electrostatic field in strongly magnetized neutron stars. The present investigation 
is mainly associated with the emission of high energy electrons from the dense metallic iron crystal 
present at the crustal region of strongly magnetized neutron stars. The study of the formation of 
plasma in a pulsar magnetosphere is a quite old but still an unresolved astrophysical issue. In the 
formation of magneto-spheric plasma, it is generally assumed that there must be an initial high energy 
electron flux from the magnetized neutron stars. At the poles of a neutron star the emitted charged 
particles flow only along the magnetic field lines. Further a rotating magnetized neutron star generates 
extremely high electro-static potential difference at the poles. The flow of high energy electrons along 
the direction of magnetic field lines and their penetration through the light cylinder is pictured with 
the current carrying conductors. Naturally, if the conductor is broken near the pulsar surface the entire 
potential difference will be developed across a thin gap, called polar gap. This is of course based on the 
assumption that above a critical height from the polar gap, because of high electrical conductivity of 
the plasma, the electric field En, parallel with the magnetic field near the poles is quenched. Further, 
a steady acceleration of electrons originating at the polar region of neutron stars, travelling along 
the field lines, will produce magnetically convertible curvature 7-rays. If these curvature 7-ray photons 
have energies > 2m e c 2 (with m e the electron rest mass and c the velocity of light) , then pairs of e~ — e + 
will be produced in enormous amount with very high efficiency near the polar gap. These produced 
e~ — e+ pairs form what is known as the magneto-spheric plasma [8ll9l fT0l[TTl[T2l[T3l[T4] . The emission 
of electrons from the polar region of neutron stars is mainly dominated by the cold emission or the 
field emission |15j . driven by electro-static force at the poles, produced by the strong magnetic field of 
rotating neutron stars. The work function plays a major role in the emission processes. Moreover, it is 
interesting to note that electrons are only emitted from the polar region of magnetized neutron stars, 
but never from some region far from the poles. A large number of theoretical and numerical techniques 
were developed to obtain analytical expression and also the numerical values of work functions of 
various materials for which the experimental values are also known. Self-consistent jellium-background 
model, embedding atom-jellium model, density functional theory (DFT) etc., were used to obtain 
numerical values for work functions [TCTTTirTg] . In a recent article the work function has been identified 
with the exchange energy of electrons within the material [Ti?Il2T)Il2T| . 

The process of extracting electrons from the outer crust region of strongly magnetized neutron stars, 
including the most exotic stellar objects, the magnetars, it requires a more or less exact description 
of the structure of matter in that region. From the knowledge of structural deformation of atoms in 
strong magnetic field; we expect that the departure from spherical nature to a cigar shape allows us 
to assume that because of high density the electron distribution around the iron nuclei at the outer 
crust region may be replaced by Wigner-Seitz type cells of approximately cylindrical in structure [221 
23,24]. We further assume that the electron gas inside the cells are strongly degenerate and at zero 
temperature (In presence of strong quantizing magnetic field the atoms get deformed from their usual 
spherical shape and also get contracted, then it is quite natural that the electron density becomes so 
high inside the cells that their chemical potential [i e 3> T, the temperature of the degenerate electron 
gas. Then the temperature of the system may be assumed to be very close to zero). It is well known 
that the presence of extraordinarily large magnetic field not only distorts the crystalline structure of 
dense metallic iron, also significantly modifies the electrical properties of such matter. As for example, 
the electrical conductivity, which is otherwise isotropic, becomes highly anisotropic in presence of 
strong quantizing magnetic field [2"51l26ll2"T] . In presence of strong magnetic field, iron crystal is highly 



Anisotropic Nature of Work Function in Strong Quantizing Magnetic Field 



3 



conducting in the direction parallel to the magnetic field, whereas flow of current in the perpendicular 
direction is severely inhibited. 

The aim of this article (i) is to show that the work function associated with the emission of electrons 
along the direction of magnetic field at the polar region of strongly magnetized neutron stars increases 
with the strength of magnetic field and (ii) it will be shown that the same quantity associated with 
the emission of electrons in the direction transverse to the magnetic field direction is infinitely large. 
Therefore the main result of this article is to show that in presence of a strong quantizing magnetic 
field, the work function becomes anisotropic (see also [201I2T] and [28 1 129 ] h We have noticed that the 
result does not depend on the dimension of the crustal matter. The scenario is very much analogous 
to the charge transport mechanism within the metal in presence of strong magnetic field. The electron 
emission process may be assumed to be a kind of charge transport process from inside the metal to 
outside. We shall show in this article that analogous to the internal charge transport, this is also 
anisotropic in presence of strong magnetic field and if the magnetic field strength is high enough, 
the emission process in the transverse direction is totally forbidden. To the best of our knowledge, 
using this simple idea for cold cathode emission [T] , the study of anisotropic nature of work function 
in presence of strong quantizing magnetic field, relevant for strongly magnetized neutron star crustal 
region, has not been studied earlier. Of course, a completely different approach, popularly known as 
density functional theory (DFT) has been used by Lai et. al. to study properties of dense neutron 
star crustal matter. They have also obtained the work function as a function of magnetic field [2"81 
29,30] (see also (31)). Surprisingly, our simple minded model gives almost same kind of magnetic field 
dependence for the work function as obtained by Lai et. al. 

We have organized this article in the following manner: in section 2, we shall obtain an approxi- 
mate analytical expression for work function associated with the emission of electrons along the field 
direction. In section 3, we shall obtain the same quantity corresponding to the emission transverse to 
the direction of magnetic field, and show that this particular component is infinitely large. Whereas in 
the last section we shall give conclusions and future prospects of this work. 



2 Longitudinal Emissions 

In this section, for the sake of completeness, we shall repeat some of the derivations presented in pQ. 
Based on this particular work we shall develop a formalism to obtain the work function associated with 
the emission of electrons along the direction of magnetic field from dense crystalline structure of mainly 
metallic iron in the outer crust region of neutron stars. It should be noted that the physical picture of 
our model is completely different from fl,. In our model calculation the arrangement of cylindrically 
deformed WS cells in the crustal region are such that the field lines and the axes of the cylinders are 
parallel with each other. To obtain the work function for both axial and transverse emission processes, 
let us consider fig.(l), where we have considered a cylindrically deformed WS cell of the crustal matter. 
Here the magnetic field is along z-axis and r-axis is orthogonal to the direction of magnetic field lines, 
and we assume azimuthal symmetric geometrical structure. To get the longitudinal part of the work 
function, we assume that an electron has come out from the cylindrically deformed cell through one of 
the plane faces, and is at the point P. The electro-static potential at Q produced by this electron is 
given by 

where AP = h, PD = r and CQ = z. This equation can also be expressed in the following form 

<tf*\r,z) = e / J (£r) exp[ T (z - h)Z]d£ 

J£=0 

for z > or < h (2) 
For the charge distribution within the distorted WS cell, the Poisson's equation is given by 



V 2 = 47ren e 



(3) 
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where <j) is the electro-static potential and n e is the electron density, given by 

n e = ^ V f^(2-6 v0 )(» 2 -m 2 -2veB) 1 / 2 (4) 

where v is the Landau quantum number with v max the upper limit (it is finite for T = 0). The factor 
2 — 5 v q within the sum takes care of singly degenerate zeroth Landau level and doubly degenerate other 
levels with v ^ and fi is the chemical potential for the electrons, given by /i = Ep — e<f){r) = constant, 
this is the so called Thomas-Fermi condition, Ep is the Fermi energy for the electrons, given by 
Ep = (p F + ml) 1 / 2 (throughout this article we have chosen h = c = 1). Hence we can re-write the 
Poisson's equation in the following form 

o 2 D vvaax 

V 2 = — Y (2 - ) [(/i + ecp) 2 -m\- 2veB\^ 2 (5) 

!/=0 

To get an analytical solution for <j> in cylindrical coordinate system (r, z), with azimuthal symmetry, 
we put, m e = 0, i.e., the kinetic energy is assumed to be high enough compared to the electron rest mass 
and v max — 0. The later is a valid approximation, provided the magnetic field is too high (> 10 14 G), so 
that the electrons occupy only their zeroth Landau level. We expect, that the above assumptions are 
not so drastic for the type of physical system considered here. Under such approximation, the Poisson's 
equation reduces to 

V 2 0«^-^ + e0) (6) 

7T 

Now defining ip = fi + e<p the Poisson's equation becomes 

9p 3 B 

w —i; = A 2 V (7) 
In cylindrical coordinate (r, 0, z), the above equation can be written as 

dr 2 r dr d 2 z 
The solution in the cylindrical coordinate system is given by 

/>oo 

1>(r, z) = / Jo(£r)o(£) exp[(^ 2 + \ 2 ) 1/2 z]d£ 
J£=o 

with z < (9) 

Here A 2 = 2e 3 _B/7r (in this connection, one should not that the parameter A used in [T] has completely 
different physical meaning) and a(£) is some unknown spectral function. Then following [I], we assume 
that there exist a fictitious secondary field in vacuum. The work function is defined as the work done 
by this field in pulling out an electron from just inside the material surface to infinity. The secondary 
field, as introduced in [T], is expressed in coherent with </>( p )(r, z) and (f>(r, z), and is given by 

4> (s) (r,z) = e / J„(£r)/(0 exp{-£z)d£ 

J£=0 

with z > (10) 

where /(£) is again some unknown spectral function. To obtain /(£), we follow [TJ and use the conti- 
nuity conditions for tangential and transverse components of electric field and the displacement vector 
respectively, given by 

E (t) = E ( P )t + E (s)t and D X _ D (p)X + d (b)X (11) 

where for the case of electron emission along magnetic field direction, the tangential part of electric 
field components are given by r derivatives of the corresponding potentials with a negative sign as the 
multiplicative factor, e.g., E^ — —d<j>/dr and similarly for others. The axial part of the components for 
displacement vectors are obtained by taking z derivatives of the corresponding potentials multiplied by 
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— 1 and the dielectric constant of the medium, which are assumed to be unity, i.e., D = —Kdcf>/dz = 
—d(f>/dz = E . Using the z derivatives for exponential functions and the relation dJo(£r)/dr — 
£Jq(£?") = —£Ji(£r) for the r derivatives for Bessel functions, we finally have following [T] 



/(0 



Then we have the secondary potential 



<j> {s) (r,z) = e Jo(Cr) 

J£=0 



1/2 



( 1 + f) 



1/2 



exp(-fh) 



(12) 



1-1 



1/2 



A 2 



1/2 



exp[-£(/i + z)]c2£ for z > 



(13) 



Hence the secondary field along axial direction at r = is given by 



E^ = -- 



s) 



1-1 



1/2 



1 



1 + | 



1/2 



expl-tth + z)m 



(14) 
(15) 



Then the force acting on an electron at z = h, which is at the verge of emission, is given by 

1/2- 



Fj"> = eE^ = e 2 



5=o 



1 



1 + & 



1 + ^ 



1/2 



exp(-2/i£)£d£ 



(16) 



The work done in pulling out an electron from just inside the metal surface to infinity is then given by 



W f = - f 
Jo 



F^dh 
1 - ( I 



1/2 



(17) 
(18) 



Substituting £ = A sinh u, we have 

e 2 

^ = y A 



exp(— 2 it) coshu du 



o 



3 e 



2e 3 B 



1/2 



Expressing the magnetic field in terms of the critical field strength, we get 



A 



1 



3 \ttB, 



2B 



(e) 



1/2 



(19) 



(20) 



where w 4.43 x 10 13 G, the typical field strength for electrons to populate their Landau levels 
in the relativistic scenario. This equation gives the variation of work function with the strength of 
magnetic field (~ f? 1 / 2 ) associated with the emission of electrons along the direction of magnetic field. 
From this expression it is also obvious that for a given magnetic field strength, if m e is replaced by m p , 
the proton mass or to/, mass of the ions, then, since to/ 3> m p m e , the work functions associated 
with their emissions along z-direction can also be expressed by the same kind of inequalities. Then 
obviously, it needs high temperature for thermo- ionic emission of protons or ions, high electric field for 
their field emissions and high frequency incident photons are essential for the photo-emission of these 
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heavier components. Since the work function oc B 1 / 2 , the emission of electrons should decrease with 
the increase in magnetic field strength. Physically, it means that the electrons become more strongly 
bound within the metals. In fig. (2) we have shown the variation of Wf, the work function associated 
with the emission of electrons along the magnetic field direction, with B, the strength of magnetic field. 
This graph clearly shows that the work function is an increasing function of magnetic field strength 
and for low field values, Wf is a few ev in magnitude, which is of the same orders of magnitudes 
with the experimentally known values for laboratory metals. Since in this approximate calculation, the 
magnetic field is assumed to be extremely high, we are therefore unable to extrapolate this model to 
very low field values. Also, in this model we can not show the variation of Wf from one kind of metal 
to another. However, a comparison of our result with fig. (4) of [55] shows that in the DFT calculation 
also the variation of work function from one kind of metal to another one is not so significant. Further, 
our result is consistent with DFT calculation. 



3 Transverse Emissions 

Let us now consider the emission of electrons in the transverse direction. As shown in fig. (I), here p is 
the position of an electron, came out through the curved surface, then the electro-static potential at q 
is given by 

^ ) fr'> ^ +(r ! n , m » (21) 

where ap — r§ , cq — r and pd = z . It can also be expressed as 

<tt*\r, z)=e J K(r - r )] exp (t*0<£ 

J£=0 

for z > or z < (22) 

whereas the form of <f>(r,z) (or the modified form ip(r,z)) and (jy- s '{r,z) remain unchanged (eqns.(6)- 
(7)). Using the relation 

J (u-v)= Jk(u)J k (v) (23) 

k— — oo 

we can express 

+ 00 

Mtr)=J [e(r + r -r )] = Jo[t(r + r )-tro]= ]T J k [£(r + r )] J fc (£r ) (24) 

k— — oo 

Then we have 

+oo 

/.-«, 

Similarly we have 

Jfc(2£r ) 

and 



-\-oo ^oo 

^(r, z) = J ^ r o) / J x U r + r o)] o(0 exp[-(£ 2 + A 2 ) 1 / 2 | z |]d£ 

k=-oo J t=° 



(25) 



ct>M(r,z) = e V J k (2£r ) / J k [^(r + r )} exp(-£ | z \)d£ (26) 



(r, z) = e V J k (£r ) / J fc [C(r + r )]/«) exp(-£ | z \)d£ (27) 

fc=-oo J t=° 
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Then from the continuity conditions on the curved surface along r direction, we have 

E J ^ r °) / J '^ r + r o)i£°(o cx pH£ 2 + a2 ) 1/2 1 z \w 

+ 00 »oo 

= e V J k (2tr ) / 4[C(r + r )]£exp(-£ | z |)d£ 

k=-oo J ?=0 
+ 00 /.00 

+ e E MCro) / 4[e(r + r )]€/(Oexp(-eU|)d€ 



(28) 



Putting z = 0, r = R and redefining the spectral function a(£) = a(£)/e, we have 

/.00 +00 

/ ^ E 4[(« + ro)^]{a(OJfl(€r ) 

- Jfc(2^r )- J*Kr )/(0}=0 
Since £g?£ is arbitrary, we have 

+00 

E 4[(* + ro)£]{a(0./fl(fr)) 

A;— — 00 

- Jfc(2^ro)-J fc (ero)/(0}=0 



(29) 



(30) 



Similarly the continuity condition along z-direction on the curved surface at z = and r = R is given 

by 



+ 00 



+ 



E Mtro)Jk[(R + ro)tMO(e + A 2 ) 1/2 ^ 

— — oc- 
+ 00 

E M2Zro)Jk[Z(R + ro)]tdS 

■ =° fc=-oc 
+ 00 



(31) 



'€=0 fc=-oc 

Here a(£) is again the redefined form. As before, from the arbitrariness of £c?£, we have 



+00 



k— — oc 

Hence we can write 
where 

and similarly 
where 



2\ 1/2 







a(0=i r (i2,ro,0 + /(0 



F(R,r ,0 



E^oo^K(^ + ro)]J fc (2gro) 



"(0(1 + 72 



G(i2,r ,0 + /(0 



Etr_oo^[^ + r-o)]J fe (2^ ) 



(32) 

(33) 
(34) 

(35) 
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Which further gives 



solving for /(£), we have 



F + f 
G + f 



A 2 



1/2 



(37) 



f(0 



r ( i + £ 



1/2 



- G 



1-1 



A 2 



1/2 



(38) 



Now using the identity 

J'k{ u ) = 7}[Jk-i(u) - Jk+i(u)] and J u (v - u) = ^ J v+k (v)J k (u) 

— OO 

we have 

Then substituting /(£) in the expression for <f>^{r,z) (eqn.(27)), we have 



(39) 



(40) 



, (r,z) = e J (£r) 



1-1 



A 2 



1/2 



The corresponding electro-static field at (R, z = 0) is given by 



or 



Ji[(R-r )(] (l + J) 



exp(— £z)d£ 



1/2 _ Jo[(fl-ro)g]J 1 (flQ 



(41) 



5=o 



( 1+ f) 



1/2 



(42) 



where we have substituted F and G from eqn.(40). The work function associated with the emission of 
electrons in the transverse direction is then given by [T] 



p OO 

W f =-e / E {s) (R,0)dr o 
Jr. 



Using the relation 



we have 



Then finally we get 



J n {u)du — 1 

/* oo /*oo -I 

/ J [(R - r )£]dr = / Ji[(i? - r Q )£}dr = - 
Jr. Jr S 



W f = -e 2 



1/2 



Jo(R£) 



l-(l+$ 



1/2 



(43) 

(44) 
(45) 

(46) 



To evaluate W/, let us consider part of this integral, given by 

1/2 



1-1 



2\V 



172^ 



(47) 
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which is the first part of work function integral (eqn.(46)). To evaluate integral I, let us put £ = A sinh 9 
and then it is trivial to show that 

I= xr~[eM0) +r (-e)? dd (48) 

4 J exp(-0) 

Obviously, the numerical value of this integral is infinity. 

The other part of the work function integral (eqn.(46)) is also diverging in nature, but these two 
divergences will not cancel each other. It is also quite obvious that in the limiting cases, R — > oo, or 
R — > 0, the diverging nature of work function for both the cases indicate that it is independent of 
the transverse dimension of the cylindrical cells. It is obvious that for R — > 0, the second term on 
the numerator of eqn.(46) vanishes, whereas the rest can be integrated analytically and found to be 
diverging in the upper limit. On the other hand, for R — > oo, the second term becomes ~ tan(£i?— 7r/4) 
(asymptotically). Which itself also diverges for certain values of the argument, but again dose not cancel 
with the infinity from the other part of this integral. Therefore in this limit also the work function 
becomes infinitely large. Such diverging nature of work function is associated with electron emission 
along the transverse direction to the external magnetic field. Although we have considered here the 
extreme case of ultra-strong magnetic field, we do expect that even for low and moderate field values, 
the work function corresponding to the electron emission in the transverse direction will be several 
orders of magnitudes larger than the corresponding longitudinal values (see also [20ll21| ). 



4 Conclusion 



In conclusion, we mention that the main purpose of this article is to show that (i) in presence of strong 
quantizing magnetic field the work function becomes anisotropic, (ii) the transverse part is infinitely 
large, (iii) the longitudinal part is finite but increases with the strength of magnetic field and finally, 
(iv) low field values of work functions are more or less consistent with the tabulated values. In the 
present work, with such simple minded model, to the best of our knowledge, the anisotropic nature of 
work function in presence of strong quantizing magnetic field is predicted for the first time. Further, the 
diverging character of work function associated with the electron emission in the transverse direction 
is also obtained for the first time, and so far our knowledge is concerned, it has not been reported in 
any published work. We have also noticed that in the low magnetic field limit (within the limitation of 
this model) the numerical values of work function are of the same orders of magnitude with the known 
laboratory data. However, we are not able to compare our results with the variation from one metal to 
another. The anisotropic nature of work function is apparently coming from the deformed cylindrical 
nature of electron distribution around iron nuclei caused by strong magnetic field at the outer crust 
region of a strongly magnetized neutron star. Now because of Ohmic decay, when the field strength 
will become low enough (< fic )j we do expect that the spherically symmetric nature of the electron 
distribution around the nuclei will be restored. As a consequence, the electron emission processes will 
no longer be affected by strong magnetic field through work function of the metallic outer crust. If the 
diverging nature of work function associated with the emission of electrons through the curved faces 
of the so called metallic atoms (they are in reality the electron distribution around the iron nuclei), 
is solely because of cylindrical deformation, which in the present model caused by the presence of 
ultra-strong magnetic field, we therefore expect that for the emission of electrons through the curved 
surface of any object having cylindrical structure, the transverse component of work function will be 
much larger compared to the longitudinal part. As a consequence the electron emission current along 
the axial direction will be large enough compared to the transverse emission current [20] . The work 
function associated with the emission of electrons along the direction of magnetic field may be used to 
obtain electron field emission current at the polar region of strongly magnetized neutron stars, which 
will help us to understand the formation of magneto-spheric plasma. 
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Fig. 2 Variation of work function (Wf) with the magnetic field strength. 



